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■SUMMARY 


After buckling, statically indeterminate trusses, beams, and other 
"strictly symmetric" structures may collapse under loads which reach limiting 
magnitudes. The current paper discusses optimal design fot prescribed values 
of these collapse loads. 


INTRODUCTION 


The principles and techniques of optimally designing struetural elements 
against buckling have been widely investigated. For example, there exists an 
extensive literature on the problem of finding the least weight design for a 
column of prescribed Euler buckling strength (see, for example, ref. 1,2,3), 
and two recent publications (ref. 4,5) deal with the analogous problem of find- 
ing the lightest beam to resist lateral buckling under prescribed loads. The 
common feature of these problems is the fact that the structures considered are 
statically determinate in the sense that the prebuckling stresses themselves 
are independent of the design. 

If the structure is indeterminate, and if the prebuckling stresses them- 
selves are therefore affected by design changes, the problem becomes vastly 
more complicated and no general optimality principles appear to have been 
developed. On the other hand, it is likely that in cases of this type the 
buckling load i.tself does not represent an important design criterion. Some 
structures buckle under decreasing loads and are therefore imperfection- 
sensitive. Others may buckle under increasing loads, and their actual strength 
is again governed by factors other than the critical buckling load. 

It has been shown that certain "strictly symmetric" types of structures 
necessarily buckle under increasing loads, and that these loads often approach 
limiting values as buckling deformations increase indefinitely. Examples of 
struetures of this kind are statically indeterminate trusses (ref. 6) or beams 
buckling laterally (ref. 7), and recent numerical (ref. 8,9) and experimental 
(ref. 10) results have confirmed the general theory (ref. 11). It may there- 
fore be realistie to study the optimal design of such structures as their 
collapse strength, rather than their buckling strength, is prescribed. The 
object of this paper is to introduce a general discussion of this problem and 
to indicate a method of solution. 
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POSTBUCKLING MODEL 


The postbuckling behavior of strictly symmetric structures has been 

described in total generality in reference 11. It can easily be visualized 

by means of a simple model consisting of a pin- jointed truss of n (say, 

n = 2) degrees of indeterminacy. If the external loads are increased by 

increasing a common load parameter X, then the "critical" load value is 

reached when the compressive force in one of the bars (say, bar 1) reaches 

the Euler value for that bar. Nevertheless, the load-carrying capacity of 

the truss is obviously not yet exhausted* While member 1 buckles under 

sensibly constant compressive force, X^continues to increase until member 2 

similarly starts to buckle. Collapse occurs when member 3 also buckles, and 

X = X then remains constant . 
c 

This simple process can be visualized within a format that is applicable 
to all strictly symmetric structures. Let S, the vector of all bar forces, 
be of the form 



+ ot S 


r ~r 


( 1 ) 


in which, for simplicity, the self-equilibrated bar force systems are 

selected so as to satisfy the orthonormality condition 
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a. 

r si 

A . E . 

1 1 



• S 

~s 



(r,s = 1,2) 


( 2 ) 


where the summation extends over all the bars and E^ represent, 

respectively, the length. Cross-sectional area, and Young's modulus of the 
i'th bar. Moreover, if Sg is the actual force system in the unbuckled 
structure, (aj. = 0), then 

S • S = 0 (r = 1,2) . (3) 

~o ~r 


In the absence of any limitations on the tensile strength of any member, 
the condition of "statical admissibility" is given by 

^ . (N^ > 0 = Euler force) , (4) 

which, in view of equation (1), becomes 

a > - N^-XS^ (i = l,2,...,n) (5) 

r r - o 

For given value of X equations (5) define a statically admissible region in 
the Uy space, whose convex boundary consists of hyperplanes whose normal 
vectors are proportional to Sj (fig. 1). The region so defined need not be 
closed. For definiteness we assume X >0 and sj < 0 (i = 1,2,3, — ip^n); in 
that case the region "shrinks" for increasing values of X. 
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For the sake of brevity we rule out the possibility of multiple buckling 
modes; then the critical value A = is reached when 


^1^0 " ’ ^1^0 ^ Ci = 2,3,...,n) . (6) 

As bar 1 buckles under constant compressive Euler force the first (i = 1) of 
equations (5) , in view of equation (6), becomes 


a = - (A - AJS^ . 
r r '• 1^ 0 

At the same time the changes in the bar chord lengths are given by 
^ A ^ " ‘^1 

^ -i ~ p ~ 2 , 3 j • • ■ > n) 


(73 


( 8 ) 


in which 6^ >0 represents the nonlinpar effect of the curvature. Hence 


I 6. = S • S - S\ = 0 (r = 1,2) , 
L r 1 ~r ~ r 1 ^ > j > 


(9) 


or, with equations (1), (2), and (3), 

a = sVcS’ (r = l,2) . 

Finally, when equation (10) is substituted into equation (7), 



( 10 ) 


(113 


confirming, once again, that strictly symmetric structures have stable points 
of bifurcatioh. ■ ; 

For A < Aj^ the origin 0 of the coordinate system in figure 1 is in 
the statically admissible region and therefore represents the actual stress 
point. At bifurcation (A = Ai3 the hyperplane passes through the origin 

and, for increasing values of A, the origin moves outside of the statically 
admissible region, while the stress point P moves with Bi. According to 
equation (103 the vector OP is parallel to the normal to B^ and, because 
of the convexity of the stable region, P is therefore closer to 0 than any 
other statically admissible point. 
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After bar 2 also buckles, point P lies on the intersection of two 
hyperplanes, and 

a = 6' + 6i . (12) 

r r 1 r 2 ^ ^ 

Finally, collapse is reached, for X = Xc, when the statically admissible 
region has shrunk to the point Pg representing the intersection of three 
hyperplanes. In that case the constant values of a^. are given by 

Cr = l,2) , (13) 

and as collapse proceeds according to 


61^ = 0)5? (oj-^») , (14) 

the collapse mechanism satisfies 

sj5^ + 8^5^ + 8^6^ = 0 (r:.l,2) . (15) 

We also note that, in general, this mode as well as the value of is 

independent of initial imperfections. 


OPTIMALITY 


For the more general case we may identify the major state of stress by 
means of 


0 = Xo + 


a a 
r-r 


The equations of compatibility are given by 



(r = 1,2, . . .,n) 


(16) 


(17) 


in which G is the compliance density with respect to a, %2 
quadratic contribution to the major Strain associated with the buckling mode 
V, and the notation implies an integral or a summation over the entire 
structure. 

The condition of equilibrium is given in variational form by 

k'^(v) Kk(5v) - g'^ £^j (v5v) = 0 , (18) 
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where k is the linear buckling strain tensor and K the stiffness density 
with respect to k. We note that both K and C are, in general, functions 
of the design variable h. 

For optimality we vary the design by replacing h by h + ^, subject to 
the condition of constant volume 

• dA • 

V = 3i^h=0. (19) 


Since the load is prescribed it follows that X = 0; nevertheless, the major 
stress system (identified by Oj.) and the buckling mode v may change. 
Variation of equations (17) and (18) then leads to 

T fdC . . . n 

~r Ldh - ?s®s ■ ~ll*-- y^J "" ^ (r - 1,2, . . . ,n) (20) 

k'^(v) K k(6v) - (v ^y) = “s'^s^ll 

t T 2 dAl • 

r(v)g^k(6v) - A^ 2A Jh (21) 

2 

A has been introduced as Lagrangian multiplier to account for 
(19). Equation (18) represents. a homogeneous eigenvalue problem, and 
(21) has therefore no solution unless the condition of integrability 

• T r T dK 2 dA T • 

“s ?s ^2^^^ ■ b ^ dh J ^ ° 

is satisfied. We note that equations (20) and (22) are similar to the equa- 
tions derived for the initial buckling case in reference 4, except for the 
last term in equation (20) representing the contribution of the postbuckling 
condition. 

Letting once again 

V = tov^ A = uA^ (w-^w) (23) 

and assuming collapse under finite load and stress conditions we obtain 

?r^2^^c^ " ° l>2,...,n) (24) 

= 0 C2S) 


in which 
equation 
equation 


T 

k ( V ) . . i\ I V 

~ ^-c"^ dh ~^-c 


dK 2 dA 

k(v^) = h — 


dh 


(26) 
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of which the first two equations represent the collapse condition, and the 
last constitutes the condition of optimality. It is noted that once again 
this optimality condition requires constant strain energy density in the 
design fibers. It is also noted that for collapse (in eontrast to initial 
buckling) the direct effect of a design change on the collapse mode via the 
compatibility conditions has disappeared". In other words, we see once again 
a parallel behavior pattern between collapse through buckling and collapse 
through perfect plasticity. 


EXAMPLE 


As an example to illustrate the theory, we consider a beam of length i 
which is fixed in its own major plane at the right end and subjected to a 
bending moment X at the simply supported left end. Collapse occurs when 


a 

c 



while the equations of equilibrium (25) assume the form 


(27) 


K,u"-a$ =0 (K.,e') ' + 0 u" = 0 (0<xsil) 

1 c c c '• 2 c C c . \ ^ J 


(28) 


where u and g represent the lateral displacement and rotation, respec- 
tively, with associated bending and torsional stiffnesses Kj and K2. In 
the development of equations (28), it is assumed that u = u" = g = 0 at both 
ends and that the effect of warping can be neglected. In terms of 3 alone 
equations (28) reduce to 


CK23^) ' + 3c = 0 (0<x <£) (29) 


The collapse condition equation (24) becomes 
£ 

X u" 3 dx = 
c c 

I 

while the optimality criterion equation (26) assumes the form 



i 


i el dx - 0 


(30) 


dK /or dK- - 

If, _£ 2 .i2 

dh \K^/ ^c dh '"e 


,2 ^ 

c dh 


(O s x s £) 


(31) 


3 

For the specific case of a thin rectangular beam, in which Kj = b h/12, 
K2=b3h/3, and A = bh, and in view of equation (29), equation (31) can be 
written in the form 
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' ( 32 ) 



which lends itself well to an iterative solution seherae. It is also 
interesting to note that equation (32) is satisfied for constant value of h 
provided 6 = sin ttx/£; this confirms the curious conclusion arrived at re- 
cently by Popelar (ref. 4) that the prismatic design represents an optimum 
for simply supported beams under constant bending moment. 


Numerical results covering equations (29) , (30) and (32) for the case 
under consideration are currently being developed. Because of the variation 
in the major bending moment, it is expected that in this case the prismatic 
beam is not optimum, and that optimal design for collapse may lead to a 
noticeable reduction in weight. 
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Figure 1,-' Redundant stress space 



